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Abstract 



A plausible explanation for the existence of additional light sterile neutrinos is 
that they correspond to modulini, fermionic partners of moduli, which propa- 
gate in new large dimensions. We discuss the phenomenological implications 
of such states and show that solar neutrino oscillation is well described by 
small angle MSW oscillation to the tower of Kaluza Klein states associated 
with the modulini. In the optimal case the recoil electron energy spectrum 
agrees precisely with the measured one, in contrast to the single sterile neu- 
trino case which is disfavoured. We also consider how all oscillation phenomena 
can be explained in a model including bulk neutrino states. In particular, we 
show that a naturally maximal mixing for atmospheric neutrinos can be easily 
obtained. 

There are several indications of neutrino mass from experiments sensitive to neutrino oscil- 
lation. Recent reports by the Super-Kamiokande collaboration |jl|, indicate that the num- 
ber of in the atmosphere is reduced, due to neutrino oscillations. These reports seem to 
be supported by the recent findings of other experiments Q, as well as by previous observa- 



tions II . The data is consistent with Ami vr ~ to 10"^) eV^, sin^ 26*^^ > 0.8 while 



Am\^^ ~ (0.5- 1.1)- 10-1° eV^, sin^ 2(9„e > 0.67, where a is /i or r. Recent measurements 



^Best fit regions for solutions to the solar neutrino deficit have been identified in 
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of the day night asymmetry and the recoil electron energy spectrum disfavour the small an- 
gle MSW solution |^ . The collaboration using the Liquid Scintillator Neutrino Detector at 
Los Alamos (LSND) has reported evidence for the appearance of v^—Vf, Q and Vy,—Vf. oscil- 
lations [^]. Interpretation of the LSND data favours the choice 0.2 eV^ < Am^ < 10 eV^, 
0.002 < sin^ IQ < 0.03. Note that if neutrinos provide a significant hot dark matter com- 
ponent, then the heavier neutrino(s) should have mass in the range ~ (1-6) eV, where the 
precise value depends on the number of neutrinos that have masses of this order of mag- 
nitude [11|. Of course, this requirement is not as acute, since there are many alternative 
ways to reproduce the observed scaling of the density fluctuations in the universe. 

If all these indications should prove to be correct there must be more than 3 light 
neutrino species because explanation of the data requires 3 different mass differences. 
This immediately raises the question why these states should be anomalously light. In the 
case of the three doublet neutrinos the underlying SU(2) gauge symmetry guarantees they 
should remain massless in the Standard Model. However additional sterile neutrinos do 
not transform under the Standard Model gauge group and thus some additional ingredient 
is needed to explain their lightness. A plausible explanation has been given recently in 
the context of theories with large additional space dimensions. 

These are models where the Standard Model (SM) fields are confined to three-branes 
which are localized in a higher-dimensional gravitation-only bulk. Bulk fermions are, 
therefore, SM singlets and are associated to Kaluza-Klein towers of 4-dimensional sterile 
neutrinos. In the context of string- or M-theory derived supersymmetric models the su- 
perpartners of moduli fields are particularly interesting candidates for light bulk fermions. 
Perturbatively, these moduli fields as well as their fermionic partners are exactly mass- 
less. Masses for the moduli and modulinos are then generated by non-perturbative effects 
which, at the same time, may also break supersymmetry. This may account for the light- 
ness of the sterile states. As an example relevant to the present paper, one can consider 
models with a string scale of order TeV and supersymmetry breaking on the observable 
brane. These supersymmetry breaking effects would be communicated to the bulk with 
a suppression of one power of the low-energy Planck scale resulting in a scale of order 
meV. Therefore, in such models, it is reasonable to expect light bulk masses for moduli 
and modulini of the order meV. 

From the experimental point of view, it seems unlikely that a single sterile neutrino is 
involved in atmospheric neutrino oscillations Moreover, as we discuss below, large mix- 
ing angles between doublet and singlet neutrinos are ruled out by the constraints following 
from the requirement that a supernova should not be rapidly cooled by sterile neutrino 
emission. The possibility that (small angle MSW) oscillations into a single sterile neutrino 
solve the solar neutrino problem also seems to be disfavoured by recent SuperKamiokande 
results The SK preliminary analysis also disfavours the standard small mixing angle 
solution in the active neutrino case, so that only the large mixing angle solution would 
be left. We will see, however, that oscillations of solar neutrinos into towers of sterile 
neutrinos are not disfavoured. On the contrary, they fit very well the presently available 
data. Furthermore the supernova bounds do allow the small mixing angle needed for this 
mechanism to work. Therefore, the brane-world scenario, besides being a theoretically 
well- motivated framework, is an experimentally viable sterile neutrino scenario, given the 
present experimental status. At the same time, it offers a small mixing alternative to the 
large mixing angle solution. 
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The connection between large extra-dimensions and neutrino physics has been dis- 
cussed in |12] and in |12, 14, 15|. In this letter we present an explicit model showing 



that the presence of a Dirac bulk mass terms can significantly affect the phenomenology 
of brane-bulk neutrino models and can open up interesting possibilities from the model- 



building point of view |15|. A general analysis of the impact of different possible mass 



terms will be given elsewhere [16|. As for solar neutrino oscillations, we will descibe new 



schemes that include a modified version of the "massless" model in |13]. A detailed fit in- 



cluding the latest SK (1117 days) and GNO data will show the beautiful agreement with 
experiment. We will then see how a maximal angle for Vr atmospheric neutrino 

oscillations can be naturally generated. The model has also room for V(, ^ oscillations 
with a Am^ larger than the atmospheric one, that can be used to accommodate the LSND 
signal. 

We first briefly describe the model from a 5-dimensional point of view. We will then 
translate it into a 4-dimensional language and study its phenomenology. We consider 
an effective 5-dimensional model with an approximate "lepton number" U(l) symmetry. 
Under this symmetry the SM leptons as well as a set of bulk spinors ^'j, where I = 
1, . . . , A^, carry charge L = iQ A 6-dimensional explicit example with a similar structure 
has been described in |15|. In this paper, we consider a 5-dimensional model with the free 
bulk action given by 

5buik = j d^xdy {^n^idA^^i - fii^i^i) . (1) 

Here, x = {x^, . . . ,x^) are the coordinates longitudinal to the brane, y = x^ is the trans- 
verse coordinate, 7"^, where ^ = . . . 4, are the 5-dimensional gamma matrices and ^/ > 
are Dirac masses. This action for the bulk fermions is the most general free one compatible 
with 5-dimensional Lorentz- and U(l)-symmetry up to a fleld redefinition. Mass terms 



breaking the Lorentz symmetry can also arise [16| but will not be considered here. 

The most naive brane-bulk coupling is that between the standard model invariants 
LiH and the five-dimensional fermion fields (or their conjugates), evaluated at the 
brane which we locate at y = 0. This coupling takes the form 

5brane = " / d^x (^^nUH + -^WnLiH) + h.c. . (2) 
Jy=o \VM5 VMs J 

Here, Li, where i = e,fi,T, are the three lepton doublets in a basis in the family space 
that diagonalizes the charged lepton mass matrix. Furthermore, H is the relevant Higgs 
field, is the charged-conjugated fermion and M5 is the effective 5-dimensional Planck 
scale. Note that the coupling h respects the U(l) symmetry. We have also included the 
coupling h^, which explicitly breaks the lepton number symmetry of the five-dimensional 
Dirac kinetic term. We will initially consider a U(l) symmetric model with = and 
later switch on this coupling as a small source of U(l) breaking^ 

^Since ^ and '^'^ have the same transformation properties under the 5-dimensional Lorentz symmetry, 
one can choose to name "i/ the bulk spinors whose lepton number is 1. Bulk fermions with a lepton number 
different from ±1 are decoupled from the SM neutrinos in the U(l) symmetric limit. 

^An alternative source of U(l) breaking is represented by 5-dimensional Majorana mass terms for the 
fermion fields. This possibility is discussed in |lq]. 
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We now describe the model from a 4-dimensional point of view. We compactify the 
total action 5brane + "^buik on a circle of radius R. We denote by ^ and r/ the Weyl 
components of the Dirac field that is = {^,1]'^), and we expand S^, rj in Kaluza- 
Klein modes as follows: 

Ci{x,y) = ^=Y^^^j{x)e-'^y/^ (2a) 
r?Kx, y) = Yl Vni{x)e'''^/^ . (2b) 

We then get the 4-dimensional mass Lagrangian 

-Cm = ^ £,nl (fJ'I - ij^j Vnl + mRUi , (3) 



where we have defined 



hiiv , ^ h%v M, 

mii = : = hiiV , mr, = , : = hjiV . (4) 

V27ri?M5 Mpi ' V2vri?M5 Mpi ^ ^ 

In eqs. (^§), f is the Higgs VEV and I'i, i = e, fj,, r, are the SM neutrino flavour eigenstates. 
We have also used the relation 2'irRM^ = (Mpi/Mg)^, which defines the effective scale M5. 
Due to the M5 suppression, the brane-bulk mixing mn is of order hji{Ms/ TeV) 0.5 ■ 
10-5 eV. 

From a 4-dimensional point of view, the model therefore involves the three SM neu- 
trinos fj, where i = e,fj,,T, and two sterile or "bulk" neutrinos for each mode 
number n G Z and each bulk fermion / = 1, . . . , A^. The two sterile modes and ry„/ 
combine into a Dirac neutrino with mass 



Mnl = ]j f^j + , (5) 

whose component ry„/ mixes with the SM neutrinos. In what follows, we will consider the 
regime defined by 

1^ COth(7ri?//7) = 2^ 2 , ^2/„2 « 1 ' (6) 

1=1 n,i,i + n /R 

in which the brane-bulk mixing can be treated perturbatively. 

In order to study the phenomenology of the model, we write the SM flavour eigenstates 
as a superposition of mass eigenstates diagonalizing the mass Lagrangian in eq. (P). In 
the perturbative limit (P), we get |15] 

= U,kh + ^Noi + V2 V ^iV„, , (7) 

n>l 

where the mass eigenstates are denoted by z)^, fc = 1, 2, 3 and Nni, n > 0. The state Nni 
has a mass M„/ given by eq. (^) and is mainly a superposition of the two bulk modes 
rini and rj^nj. It turns out to be the left-handed component of a Dirac spinor whose 
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right-handed component is mainly a superposition of the fields and £,-ni- The three 
Majorana neutrinos where k = 1,2,3, are exactly massless in the limit of an exact 
U(l) symmetry. They acquire mass once the lepton number violating couplings h'j- are 
switched on^. Their masses and mixings Uik with the SM neutrinos can be obtained by 
diagonalizing the Majorana mass matrix m', approximately given by 

(miimjj + mj-mijj fij 
m^j = - 2^ -^^ 2 I 2/P2 ^ {mnmjj + m%mij) coth(7ri?/u/) . (8) 

Note that this mass matrix has a "see-saw" structure, where the role of the heavy mass 
is being played by (7ri?coth(7ri?/i/))~^, that is, roughly speaking, by the lightest between 
/i/ and {ttR)~^. 

Eqs. determine the phenomenology of the model and express the phenomeno- 

logically relevant quantities in terms of its parameters. The bulk masses fij coincide with 
the masses of the lightest bulk states Nqj while the splitting of the heavier bulk modes fn/ 
is determined by the inverse radius 1/R. Once the mass M„/ of the mode Nni is known, its 
mixing with the SM neutrino Vi is determined by the brane-bulk mass term m/j. Finally, 
with these parameters being fixed, the masses and mixings of the light neutrinos are 
determined by the lepton number violating parameters m^. 

Notice the difference from models where the bulk mass terms /i/ are absent. In such 
models, or in cases where <C m, the massless neutrinos decouple from the SM neutrino 
oscillations and the mass of the lightest neutrinos is determined by the same parameters 
(namely the brane-bulk masses m) responsible for their mixing with the flavour eigenstates. 
In contrast, if /i > m, the light states give a significant contribution to neutrino oscillations. 
This effect has been widely neglected in the literature so far. In particular, it offers 
the interesting possibility of associating some of the oscillation signals with oscillations 
into light neutrinos and some with oscillations into bulk neutrinos. In the perturbative 
regime we are interested in, the oscillations of SM neutrinos mainly involve the three light 
Majorana states, Vi ~ ^^^=1 Uik^k, with the approximately unitary matrix U playing the 
role of the MNS matrix. On the other hand, the oscillations into sterile neutrinos are 
associated with the small component m*j^/MojNQj + ^/2'^^y-^^m*j-/MniNni- Despite its 
smallness, this bulk component can have the important role of depleting the solar neutrino 
flux through a small mixing angle MSW effect. Naturally maximal v^- oscillations 

can then be easily obtained from the light mixing matrix (^) . The bulk physics is therefore 
(indirectly) involved in atmospheric neutrino oscillations as well, since it determines the 
light mass matrix through the see-saw formula (^). The LSND signal can also be easily 
accounted for by oscillations into light states. 

We now apply the results just obtained to the description of neutrino oscillation in 
realistic models. We first consider the one generation case of mixing between the electron 
neutrino Ve and the tower of states associated with a single bulk fermion of mass /ig, in the 
limit of unbroken U(l). Later on, we will see how the results obtained can be embedded 
into a three family scheme. 

*Once U(l) is broken, the mixing with the SM neutrinos and the heavy modes Nni (and the heavy 
modes themselves) also receives corrections. 
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Figure 1: Fit of total rates in the siv? 29Q-Am'j} plane for different values of 1/R. From 
right to left, the four patches correspond to 1/R = (20,5,2,0.5) • 10~^eV. The lighter 
(darker) regions represent the 99% (90%) CL. 

Figure 2: Survival probabilities for fi^ — 2.1-10^^ eV and {1/R,me) = (20, 0.058) -10^^ eV 
(dotted line), {l/R,me) = (5,0.035) • 10~3 eV (darker solid line), {1/R,me) = (2,0.023) • 
10"^ eV (dashed line) and {l/R,me) = (0.5,0.012) • 10"^ eV (lighter solid line). 



The properties of the system are completely determined by three parameters: the 
inverse radius 1/R, the mass of the bulk fermion //g and the brane-bulk mixing mass 
mg. The latter can be considered to be positive without loss of generality. The electron 
neutrino is mainly made of the light eigenstate P and has a small mixing 9q ~ rrie/ne 
with the zero mode z^o and a small mixing 0^ — \/2?7ie/il^j^ with the 7z-th mode iVj^, where 
n > 1. The squared mass difference between the mass of the light and the n-th mode is 
Ami = ^^l + n'^/R^■ 

For each value of 1/R, the agreement of the model with the total rates measured 
in the SuperKamiokande |17, |8|, Homestake [^] and Gallium |l^, 21 1 experiments 
depends on /ig and rUe, or equivalently on the zero-mode squared mass difference Attiq 
and mixing angle ^o- In Fig- regions in the sin^ lO^-ArnQ plane allowed by a fit of 
the total rates at 90% and 99% confidence level are shown for four values of 1/R, namely 
1/R= (20, 5, 2, 0.5) -10-3 eV. The BP98 solar model @ and the BPOO electron density § 
have been assumed. The figure shows that the total rates can be reproduced at 90% CL 
for each value of 1/R. The best fit is obtained for 1/i? = 5 • 10"'^ eV which, as we will see, 
also gives the best fit of the SK recoil energy spectrum. The second best fit is the one for 
1/-R = 20 • 10~^ eV, which is however problematic for the energy spectrum. The two cases 
of small 1/R give an acceptable fit for both the total rates and the energy spectrum. 

In order to understand the main features of Fig. [l|, in Fig. |2| we show the survival 
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probabilities associated to the four values of 1/R. The other two parameters //g and rrif. 
have been chosen in the corresponding best fit regions in Fig. |l]. In particular, for all 
curves we have used /Ug — 2.1 • 10~^eV. In fact, as Fig. || shows, the allowed range for 
//e is independent of \/R in first approximation. This is because /Xg defines the energy 
E ~ /ig/(2I4;) at which solar neutrinos start undergoing resonant conversion. Here, is 
the matter induced potential in the core of the sun. The parameter /ie niust therefore be 
in the range ~ (2 ~ 3) • 10~^ eV to ensure that Beryllium neutrinos are converted but 
pp- neutrinos are not. 

Let us first consider the "large case 1/R = 0.02 eV, illustrated by the dotted 

line in Fig. |2|. In this case, the resonant mixing between the electron doublet neutrino 
and the sterile tower will be significant only for the lowest state Nq, whose squared mass 
difference with the lightest state is Amg = /ig. This is because 2EVc < Ami — /^e + 1/-^^ 
for any solar neutrino energy. The phenomenology of this case is the same as the one in 
models with a single sterile neutrino — the Kaluza Klein origin makes no difference. The 
precise value of does not really matter as long as 1/R^ 7^g. On the other hand, mg 
determines the mixing angle 6o and hence the slope of the curve at the larger energies in 
Fig. |2|. The chosen value is mg ~ 0.58 • 10"^ eV and it corresponds to sin^ 26^0 ~ 0.003. 

The darker solid line represents the particularly interesting case 1/R ~ 5 • 10^^ eV 
in which, besides the light state z), three bulk states, A'^O) -^i and N2, are involved in 
solar neutrino oscillation. The survival probability falls after each resonant energy since 
one more level has to be crossed by the neutrino on his way to the earth to be detected 
as an electron neutrino. This is compensated by the stronger raise between subsequent 
resonances, so that the prediction for the measured total rates is approximately the same 
as in the previous case. Such an higher slope is obtained by using a smaller mixing, namely 
sin^ 2^0 — 0.001, corresponding to m-g ~ 0.35 • 10~^eV. It is then clear why the allowed 
regions in Fig. || move towards the left for smaller 1/R. 

The dashed line corresponds to /igi? ~ 1, rrig ~ 0.23 • 10~^eV. This situation is 



similar to that analyzed by Smirnov and Dvali |13]. However, the splitting between the 
low- lying levels is not constant as in their model but is rather governed by eq. (^). As 
a consequence, the energy dependence of the survival probability near its minimum is 
slightly different. Finally, the light solid line corresponds to a case in which fx^R > 1: 
1/R ~ 0.5 • 10"'^ eV, m-g ~ 0.12 • lO^^eV. A relatively large number of levels can now 
undergo resonant conversion. The single resonances are almost invisible because of the 
averaging over the neutrino production point in the sun. 

So far we have compared the prediction of the model with the data on total rates 
only. Recently, the SK collaboration has claimed that their data on the recoil energy 
spectrum and day-night asymmetry allow to disfavour oscillations of solar neutrinos into 
a single sterile neutrino and into an active neutrino through the small mixing angle MSW 
mechanism. This conclusion depends on the way rates and energy spectrum are combined 
in a global fit. In what follows, we will quote confidence levels for total rates and energy 
spectrum separately. Whatever is the procedure followed, we now show that oscillations 
into more than one sterile neutrino are not only still allowed, but also fit particularly well 
the measured spectrum. The energy spectra associated with the four survival probabilities 
in Fig. ^ are shown in Fig. ^. The predictions are compared with the latest data |2^, ^]|^. 



'As the data points have been graphically reduced from ref. |8|, there could be small deviations of the 
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Figure 3: Energy spectra associated to the four survival probabilities in Fig. [J 

Figure 4: Possible choice of parameters that allows to account for the solar, atmospheric 
and LSND oscillation signals. The scales corresponding to the mainly active neutrinos 
^1,2,3 CLnd to the mainly sterile neutrinos N^n, N^ri, o,re shown on the right hand side. 
Also shown are the scales of the lepton number conserving parameters me,^,T cmd of the 
lepton number violating parameters rn'^ . 



The dotted histogram again corresponds to the large \/R regime equivalent to the case of 
a single sterile neutrino. Since it has the highest slope, it gives the worst fit of the energy 
spectrum and is unacceptable at 95% CL (free normalization). However, it is possible to 
find a better agreement with the energy spectrum for mixing angles at the lower border of 
the 90% allowed region in Fig. P. The day- night asymmetry disfavours such lower angles 
for active neutrino oscillations, but does not play a crucial role in the sterile case. 

The presence of additional sterile states reduces the slope of the predicted energy spec- 
trum. This is illustrated by the (almost coincident) dashed and light solid histograms, 
corresponding to the two small 1/R cases in Fig. |2| The lower overall slope of the proba- 
bilities is also clearly visible there (the asymptotic behaviour for large E has been discussed 
in |]l^]). The integration over the neutrino energy and the convolution with the detector 
resolution smear out the difference between the survival probabilities in Fig. |2[ Quantita- 
tively, the dashed and dotted lines agree pretty well with the measured spectrum. For the 
1/R = 0.002 eV histogram, for example, the fits of both total rates and energy spectrum 



experimental points and error bars in Fig. ^ from the original ones. 

^In that region a global fit of total rates, energy spectrum and day night asymmetry still allows oscilla- 
tions into a single sterile neutrino. This is because the information about the relatively bad fit of the total 
rates is lost in points of the parameter space where the for the spectrum is accidentally lower than the 
expected value for a "true" theory. 
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are acceptable at 50% CL. 

Finally, the darker solid line in Fig. |3| shows that the corresponding survival probability 
in Fig. § beautifully fits the energy spectrum. In this case, the fit of total rates and energy 
spectrum are already acceptable at the 25% and 10% CL respectively. The total is 11 
for 21 data points. 



Unlike the energy spectrum, the day-night asymmetry |24, g] does not discriminate 
between the different possibilities discussed. In fact, the day-night asymmetry turns out 
to be small compared to the experimental errors in all cases considered above. Small earth 
effects affect the "midnight" bin only, where they can suppress the neutrino flux by a factor 
ranging from to 0.5%. The day-night asymmetry for sterile neutrinos is smaller than for 
active ones for at least two reasons. First, the matter induced potential in the earth is 
smaller than for active neutrinos by more than a factor two. As a consequence, the MSW 
resonance takes place at higher energies, where the neutrino flux is smaller. Moreover, 
a larger resonant energy corresponds to a larger resonant wavelength and therefore, for 
a given length travelled in the earth, to a smaller oscillating term. We consider a small 
day-night asymmetry in reasonably good agreement with the present data, given that SK 
finds only a 1.2o" deviation from zero asymmetry. 

The scenario discussed leads to various experimental signatures. First of all, a crucial 



test will be provided by the neutral/charged current ratio measurement in SNO |25]. 
Moreover, different sterile neutrino scenarios can be distinguished by the shape of the 
energy spectrum. As discussed, the present data already favours a flatter spectrum, which 
can be obtained if more than one sterile neutrino takes part in the oscillations. The charged 
current spectrum that will be measured by SNO could provide additional information. 
Finally, the measurement of a day-night asymmetry significantly lower than zero would 
disfavour the scenario with many sterile neutrinos. 

Before discussing how this fits into a complete three generation model we briefly discuss 
the bounds on bulk neutrinos that come from the requirement that cooling in supernova 
due to bulk sterile neutrino emission should not be greater than the energy carried off by 



ordinary neutrinos. A more detailed discussion will be given in reference [26|. The most 
efficient mechanism for producing sterile neutrinos in supernova is resonant conversion |14]. 
The neutrino effective potential, V, in the supernova core is of order 10 eV, giving a 
resonant mass Mn'^n/R:::^ {2VEYI'^ ~ 5 • 10^ eV (E^/lOO MeV)!/^ for a neutrino energy 
E. The effect of oscillation by resonant conversion can be reliably estimated for small 
mixing angles, because the width of each resonance is smaller than the separation between 
resonances. Thus the survival probability for a standard neutrino produced in the core is 
given by the product of the survival probabilities in crossing each resonance Pyy ~ n„P„, 
where P„ ~ e"'"'^"/^, 

- dV/dr ' ^ ' 



and m is the brane bulk mass. Following |]14[, we approximate the density profile by 
p{r) oc e~^/'''=°"'°. This gives V/{dV/dr) = rcore — 10^ cm and hence 7„ ^ m^/(10~^eV^) • 
(lOOMeV/S). To determine the energy loss due to coherent processes we must determine 
the survival probability of the active neutrino after production by the process e^e~ — > z^rQ. 

''Electrons and positrons are in thermal equilibrium with the neutron core through electromagnetic 
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The initial energy of the neutrinos is roughly eight times the temperature. Subsequent 
scattering of the neutrino by thermalisation processes rapidly reduce its energy. The pro- 
duction rate is very rapid and neutrinos quickly reach their thermal equilibrium abundance 
(at T = lOMeV this is achieved in 0.1ms). Thereafter thermal processes limit the neutri- 
nos to their thermal abundance. To good approximation the resonant conversion occurs 
only over the first few thermal path lengths, before annihilation processes eliminate the 
neutrino. Thus, after production, the survival probability will be 1 — (7r/2)(5n7„, where 6n 
is the number of levels crossed in one path length, 6n ~ {5p/2p)MnR, Sp/p = A(T, E)/rcore, 
where X{T,Ey) is the mean path length for a neutrino of energy E before annihilation. 
Putting all this together one finds the energy carried off by sterile neutrinos (again the 
dominant losses occur at the high initial temperature during the first second) 

l^Q'(E'^)y,rn„<lE (10) 



2 2 2 7 dE 
- 8.10e^f!^V^V''ergs- 

where V is the volume of the core and dQe{E, T)/dE is the e"'"e~ — > vV rate given 

= / - ^-"^ ■ (w)' ■ 

Requiring that e be less than the energy carried off by neutrinos gives the bound 

m^R\ o/30MeV\^-^ , , 

For the cases discussed above with {l/R,me) = (20,0.058) • lO^^gy^ (5,0.035) • lO'^ eV, 
(2,0.023) • 10-3 eV and (0.5,0.012) • lO'^eV, this bound is satisfied for T quite close to 
30MeV. Thus we conclude that a solution to the solar neutrino problem involving oscil- 
lation to bulk Kaluza Klein states is just consistent with a supernova bound^. However 
it is clear that the large angle MSW and vacuum solutions are ruled out by the supernova 
bound as is the large angle atmospheric oscillation into sterile bulk states. 



The discussion above shows that the mixing of the electron neutrino with a bulk fermion 
provides a "small mixing angle" sterile neutrino solution of the solar neutrino problem in 
good agreement with the present data. We now describe a model capable of explaining 
both solar and atmosheric oscillations in which the solar oscillation is predominantly due 
to oscillation from to a sterile KK tower as explained above and atmospheric neutrino 

processes which may be considered instantaneous by comparison with the weak processes generating sterile 
neutrinos. 

*We have allowed for a factor of 1/2 since only neutrinos undergo resonant conversion and another 
factor of 1/2 since the neutrino (which suffers many elastic collisions in a length A(r, E^)) does not always 
travel normal to the density contours 

^This conclusion differs from that of reference | |l^ . This is because we have used Sn rather than the 
total number of levels crossed in emerging from the supernova when estimating the survival probability. 
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oscillation is predominantly due to oscillation from to z^r0- Indeed, in the perturba- 
tive regime specified in eq. @, the bulk component of the SM neutrinos is too small to 
significantly affect the depletion of atmospheric neutrinos. 

The light masses and the MNS mixing matrix U are determined by the light mass 
matrix m' in (S). The simplest texture leading to maximal i^t oscillations is 



m' = \ €e^, €f, 1 m^r , (12) 















I J 
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where the e parameters are much smaller than 1. In the limit in which the e parameters 
are vanishing, the muon and tau neutrinos are superpositions of two degenerate states 
i>2 and z^s with mass |m^T-|- The small e^^r parameters (assumed real for simplicity) are 
necessary to generate a mass splitting 

Am23 - e Im^^rl, e = e^ + er 

between the degenerate states and therefore an "atmospheric" squared mass difference 
Am^rpj^ = Am^g ~ 2e|m^T-p. The corresponding mixing angle O^t is almost maximal, 
s\v?26^r — 1 — (e^ — ^t)"^/^- The ee^x, (-er parameters are constrained to be small by the 
Ve disappearance experiments. 

Notice that the simple mechanism used here to generate a maximal ^ mixing 
does not work in a three neutrino scenario aiming at explaining at the same time the 
solar neutrino data. This is because the other two squared mass differences available turn 
out to be larger than the atmospheric one, unless all neutrino masses are very nearly 
degenerate |2^. As a consequence, for non-degenerate neutrinos, there is no room for the 
small Am^ required by the solar data. 

If the texture (|l^ ) accounts for atmospheric neutrino oscillations, the two degenerate 
neutrinos can provide a cosmologically significant source of dark matter. Moreover, the 
electron neutrino can oscillate into muon or tau neutrinos with a squared mass difference 
Am^ > Am^rpj^ and small amplitudes lieeji i fier)^ respectively. It is obviously tempting 
to associate such oscillations with the LSND signal. The MiniBooNE experiment will test 
this possibility and will cover a relevant portion of the parameter space for such short- 
baseline oscillations. A short-baseline neutrino factory could further extend the sensitivity 
in the mixing parameter |25]. 

Let us now see how the texture (p^) can be obtained. We consider a model with three 
towers of sterile neutrinos to match the three families of doublet neutrinos. One tower 
is coupled to the electron neutrino and generates solar neutrino oscillation as discussed 
in the previous Section. The remaining two towers do not play a direct role in neutrino 
oscillation and must be heavy. We assume that the bulk fermions ^'e,/i,T have the same 
individual lepton numbers -^>e,/i,T as the corresponding neutrinos. The zero order form of 
eq. (^) in which the e parameters vanish then simply follows from a first stage of lepton 
number breaking which leaves — L^- and Lg unbroken. A further breaking of — Lj- 
and Le finally generates the small entries in (IT^). For simplicity we break lepton number 
through the brane-bulk couplings in eq. (pi) only. 



""Supernova bounds prohibit large mixing with sterile neutrinos so this is the only viable possibility. 
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The most general — and Lg conserving structure of the mass matrices m and 
is 



(13) 



where rue^^^r can be made positive but mj^^ and are in general complex. Moreover, 
the Dirac mass matrix of the bulk fields is diagonal, that is /i = Diag(/Xe, /x^, /i,-). The 
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light mass matrix resulting from eq. (^) is then as given in eq. (12) with vanishing e terms 
and nifj^r specified by 



m 



^r — — vri? {m^nf^r coth(7ri?^^) + mrrrir^ coth(7ri?/iT-)) . (14) 
The LSND signal, if confirmed, would determine the size of m^r as jw^Tp — ^'^lsnD' 

The 

conditions eq. @ for being in the perturbative regime become (taking also into account 
the U(l)-breaking) 

ttRIth P 

^—^coth{TTRne)<^l (15a) 

7ri?lm„P TriJlml;,,^ 

coth(7ri?^„) + ^ coth{TTRfir) < 1 (15b) 

coth(7ri?^^) + coth(7ri?^^) < 1 . (15c) 

At this stage, the electron family is decoupled from the other two. The solar neutrino 
phenomenology is therefore given by the one family model described previously, where fig 
and me are as discussed. As reference values for these parameters, we will, in the following, 
consider the best fit values fie ^ 2.1 • IQ-^eV, 1/i? = 5 • IQ-^eV, rUe 0.35 • IQ-'^eV. 
To complete the picture, we only need to break — Lr in order to generate Am^rpj^ 
and to break Lg in order to account for the LSND signal. This is easily accomplished by 
switching on the — Lr and violating entries in either m or in eq. dlD. This also 
adds to the electron neutrino a small and ^r component. Such a component, however, 
does not affect the solar neutrino discussion since the squared mass difference with the 
lowest mode is too large to give rise to resonant conversion. 

The possible combinations of parameters leading to viable models are numerous. With 
the only purpose of exhibiting an example, let us consider a case in which the lepton 
number conserving parameters me, m^ and nir are hierarchical, with a gap of two orders 
of magnitude between "near" families and me taking the reference values above. Similarly, 
we take for /ie the value required by solar neutrino oscillations. The values for fj,^ and fj,r 
do not affect the particle phenomenology of the model significantly as long as they are 
both larger than 1/R and smaller than the string scale. This is because, from eq. (|8|), 
the heavy mass scale that governs the see-saw suppression of the light masses is always 
given by 1/R for /i^ and fir in this range. As a consequence, we are free to choose for 
example and fir safely above \/2EV^ where V is the matter potential in the core of 
a supernova and E is the relevant neutrino energy as discussed above. This avoids any 
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possible problem with the supernova bound eq. (^) due to conversion of /i or r neutrinos 
into their associated Kaluza-Klein states^ 

In order to generate the LSND squared mass difference, should be given by 



where we have also assumed |mj^^| ~ Im^^l- The atmospheric squared mass difference 
and the LSND mixing angle ^lsnd can, for example, be generated by a m^^'^ 
m^g^'^ri^e mass term respectively, if 

Am^rpj^ J c sin26'LSND c 
"T-T = T^;^ — o '"T-Tii and m_„ = ^ m_,, . 

^^''''LSND 

Also, the perturbative conditions, eqs. (p^), are satisfied. 

The hierarchy of the various scales in this example is illustrated in Fig. ^ for A/zi-Lgj^j-, 
leV^ and Am^rpj^ ~ 3 • 10"'^ eV^. Notice that the brane-bulk masses are given by rrii ~ 
/ij(Ms/ TeV) 0.5 • 10~^eV, where hi, i = e,fi,T, is the brane coupling of the i family. 
Brane-bulk masses larger than 10~^ eV, as required in the above example, can arise from 
a string scale larger than ITeV or, alternatively, from large couplings hi. 

To summarise, we have studied the possibility that doublet neutrinos may oscillate 
into the sterile tower of Kaluza Klein states associated with bulk neutrinos, modulini, 
propagating in new large dimensions. The most stringent constraint comes from supernova 
and we found that large mixing angles between doublet neutrinos and light singlet bulk 
states are ruled out but that small mixing angles, of the order needed to explain solar 
neutrino oscillation, are allowed. A detailed study of the phenomenology associated with 
small angle MSW oscillation to a sterile Kaluza Klein tower shows that there is a much 
larger range of possibilities than has hitherto been explored. These arise if one allows for 
bulk masses in addition to the usual Kaluza Klein masses and allow for a range of different 
spacing between the bulk states. As a result one may go from the limit in which only a 
single bulk state is involved in the resonant conversion to the case where a continuum 
contributes. The former gives the same phenomenology as a single sterile neutrino and 
is disfavoured by the SK data for the energy spectrum and the total rates. However, if 
more bulk states contribute the deviation from the observed energy spectra reduces and an 
excellent fit to the data is possible for a level spacing in which three KK modes contribute 
significantly to the resonant oscillation. This model can readily be extended to describe all 
present indications of neutrino mass, including maximal mixing for atmospheric neutrinos 
in a natural way and the LSND signal. It can also give rise to a cosmologically significant 
dark matter component. 



^^It could even be possible that and fir are close to the string scale. As a consequence, we have a model 
in which two of the bulk states receive large masses of order TeV while one of them is small and of order 
meV. In the context of string models, the former masses could be attributed to non-perturbative effects in 
the bulk and the latter to non-perturbative effects on the brane which are gravitationally communicated 
to the bulk as discussed in the introduction. 
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